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ABSTRACT 

We report on the first detection of a pre-supernova (SN) evolution in the X-ray regime. The results 
are based on ROSAT observations of SN 1993J ranging from six days to five years after the outburst. 
The X-ray observations are used to probe the SN shell interaction with the ambient circumstellar matter 
(CSM). After exploring various scenarios that might explain the observed X-ray lightcurve with a i~ - 27 
rate of decline, we present a coherent picture in terms of the interaction of the SN shock front with the 
CSM deposited by the progenitor's stellar wind. During the observed period, the SN shell has reached a 
radius of 3 x 10 17 cm from the site of the explosion, corresponding to ~ 10 4 years in the progenitors stellar 
wind history. Our analysis shows that the mass-loss rate of the progenitor has decreased constantly from 
M = 4 x 1CU 4 to 4 x 1CU 5 M Q yr _1 (u w /T0 km s _1 ) during the late stage of the evolution. Assuming 
a spherically symmetric expansion, the circumstellar matter density profile is found to be significantly 



flatter (p c 



-1.63 



) than expected for a constant mass-loss rate and constant wind velocity profile 



(r ). The observed evolution either reflects a decrease in the mass-loss rate, an increase in the wind 
speed or a combination of both, indicating that the progenitor likely was making a transition from the 
red to the blue supergiant phase during the late stage of its evolution. 



Subject headings: supernovae: individual (SN 1993J) 
1993J) — X-rays: ISM 

1. INTRODUCTION 

The interaction of a supernova (SN) with circumstellar 
medium (CSM) produces a fast shock wave in the CSM 
and a reverse shock wave into the outer supernova ejecta. 
Two characteristic regions of X-ray emission are produced 
by the interaction: the forward shock wave in the CSM at 
T ~ 10 4 km s _1 produces gas with ~ 10 9 K, and the re- 
verse shock wave in the supernova ejecta ~ 10 3 km s _1 less 
than the forward shock produces gas with T ~ 10 7 K. The 
reverse shock is formed where the freely expanding super- 
nova ejecta catches up with the CSM shocked by the blast 
wave. At early times, the reverse shock front is radiative 
and a dense, cool (T < 10 4 K) shell can build up down- 
stream from the radiating region (Chevalier & Fransson 
1994). The dense shell can absorb X-rays from the reverse 
shock region and reprocess them to lower energies. 

This scenario was supported by observations of Type lib 
SN 1993J with ROSAT (Zimmermann et al. 1994, 1996) 
and with ASCA (Kohmura et al. 1994). The early X- 
ray spectrum at day 6 was hard, with T ~ 10 8 5 K. At 
day ~ 200, a softer component with T ~ 10 7 K domi- 
nated. The emergence of the softer spectrum could be at- 
tributed to the decreased absorption by a cool shell (Frans- 
son, Lundqvist & Chevalier 1996; hereafter: FLC96). 

2. theoretical background 

The thermal X-ray luminosity L x of the shock heated 
plasma is expressed by the product of the emission mea- 
sure, EM, and the cooling function, A(T, Z, AE), where 
T is the plasma temperature, Z represents the elemen- 
tal abundance distribution, and AE is the X-ray en- 
ergy bandwidth. For spherically symmetric conditions 

EM = 47r f^ o n 2 r 2 dr, where r, the radial coordinate, runs 
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from Rq to R, the current outer boundary of the shocked 
matter which has an electron density n c . Assuming a con- 
stant supernova shock wave speed v s , R = v s t, where t 
represents the time elapsed since the explosion. If the am- 
bient matter density p is dominated by a wind blown by 
the progenitor star of the supernova the continuity equa- 
tion requires M = 4nr 2 pv w , with M the mass loss rate 
and v w the wind speed. The X-ray luminosity can then 
be expressed as L x = ^A(T) x (M/v w ) 2 x (v s i) -1 . We 
can hence use the X-ray luminosity to measure the ratio 
M/v w of the progenitor star. 

Each X-ray measurement at t is related to the corre- 
sponding distance from the site of the explosion. This site 
had been reached by the wind at a time depending on v w , 
or the age of the wind t w = tv a /v w . Usually, v s ^> v w 
so that with t only a few years we can look back in time 
quite an appreciably large time span in the evolution of 
the progenitor's wind. Assuming that u w did not change 
over t w we can even directly measure the mass loss rate 
back in time. Whether v w = const is a reasonable assump- 
tion will be discussed below. Integration of the mass-loss 
rate along the path of the expanding shell gives the mean 
density inside a sphere of radius r. For a constant wind 
velocity v w and mass-loss rate M, a p csm — po( r / r o)~ s 
profile with s — 2 is expected. 

After the expanding shell has become optically thin, it 
is expected that emission from the SN ejecta itself, heated 
by the reverse shock, dominates the X-ray output of the 
interaction regions due to its higher emission measure and 
higher density. For a uniformly expanding ejecta the den- 
sity structure is a function of its expansion velocity, v, and 
the time after the explosion, t: p an = po(t/to)~ 3 (v/vo)~ n 
with po the ejecta density at time to and velocity vq 



2 



Immler, Aschenbach & Wang 



30 



20 



10 



H 1 — I I I I | 1 1 1 1 — I I I I | 



1 1 1 — I I I I | 1 — 

□ PSPC 
+ HRI 



w 

V) 

■*-> 
o 

CO 

I 

o 

-!-> 

tC 

S-, 

■<-> 

c 

O 

° 10 100 1000 

time since explosion (days) 

Fig. 1. — ROSAT (0.5-2 keV band) X-ray lightcurve of SN 1993 J. Boxes mark PSPC data, HRI data are indicated by crosses. Error bars 
are la statistical errors. The dotted line illustrates a t~ - 27 rate of decline. 
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(FLC96). For a red supergiant progenitor, the power-law 
is rather steep with index n ~ 20 (Shigcyama et al. 1994; 
Baron, Hauschildt & Branch 1994; Suzuki & Nomoto 1995; 
FLC96). For constant n, the radius of the discontinuity 
surface between the forward and the reverse shock evolves 
in time t with R c oc t m where m = (n — 3)/(n — s) is the 
deceleration parameter. 

3. X-RAY OBSERVATIONS AND ANALYSIS 

SN 1993J in M81 was observed in 5 individual pointings 
with the Position Sensitive Proportional Counter (PSPC) 
and in 10 individual pointings with the High Resolution 
Imager (HRI) onboard ROSAT (Triimpcr 1983). Total in- 
tegration times with the PSPC and HRI instruments are 
52.2 ks and 140.3 ks, respectively. Assuming an initial ex- 
plosion on March 28, 1993 (Ripero 1993), the observations 
cover a period of 6 to 1 821 days after the outburst. A de- 
tailed description of the data calibration and analysis can 
be found in Immler & Wang (2001). 

In order to investigate the rate of decline for SN 1993 J, 
we constructed a combined PSPC+HRI lightcurve, which 
is presented in Fig. 1. HRI data were binned into 11 
continuous observation blocks with 5-20 ks exposure time 
each to obtain satisfactory counting statistics. PSPC data 
were binned into 9 intervals with 2-18 ks integration time 
each. The time-dependent background was determined 
by normalizing the total background map of the HRI and 
PSPC images according to the total source-removed count 
rate in each exposure interval. Background subtracted 
source counts were extracted within the 90% radii around 
the fixed position of SN 1993J. In order to reduce back- 
ground due to UV emission and cosmic rays, only HRI PI 
channels 2-10 were used. 

An effective 0.5-2 keV band cooling function of A = 

3 x 10~ 23 ergs cm~ 3 s _1 is adopted for an assumed 
optically thin thermal plasma with a temperature of 
10 7 K (Raymond, Cox & Smith 1976). The equiva- 
lent count rate to (unabsorbed) flux conversion factor is 

4 x 10~ n ( ergs cm~ 2 s _1 )/(counts s _1 ) for a Galacti col- 
umn density of Ah = 4 x 10 20 cm~ 2 (Dickey & Lockman 
1990). ROSAT HRI and PSPC count rates were converted 



using a factor of 3 appropriate for the assumed source 
spectrum (cf. Immler & Wang 2001). HRI count rate 
to flux conversion factors as a function of spectral model 
parameters are given in Fig. 5 in Wang, Immler & Pictsch 
(1999). The uncertainty of the HRI (PSPC) conversion 
factors for optically thin thermal spectra with tempera- 
ture in the range 10 7 -10 8 5 K (cf. §1) is 22% (19%). An 
increase in absorbing column from 10 20 to 10 21 cm~ 2 in- 
creases the unabsorbed source fluxes by ~ 15%. Assuming 
a 0.86 keV thermal Brcmsstrahlung spectrum (correspond- 
ing to a temperature of 10 7 K) instead of a Raymond & 
Smith thermal plasma reduces the HRI (PSPC) conversion 
factor by 4% (1%). 

4. DISCUSSION 
4.1. X-ray emission from the shock-heated stellar wind 

We will first discuss the observed X-ray lightcurve in 
the context of the shock-heated stellar wind model and 
explore other scenarios later. In order to derive the stel- 
lar wind age prior to the outburst, an initial wind ve- 
locity of v w — 10 km s" 1 and shock front velocity of 
v s = 19 900 km s _1 are assumed (Bartel et al. 1994). The 
choice for the stellar wind speed is justified by the fact 
that SN 1993 J is a Type lib SN with a massive (~ 15M Q 
ZAMS) red supergiant progenitor (e.g. Podsiadlowski et 
al. 1993). Typical outflow velocities for these progenitors 
are in the range of 3 ^ v w ^ 30 km s _1 , peaking at ~ 10- 
15 km s _1 (e.g. Barnbaum, Kastner & Zuckermann 1991). 

Since the shock front catches up with the wind deposited 
by the progenitor with a speed of up to <~ 2 000 x larger, 
the interaction front can be used to probe the stellar wind 
history over a period of ~ 10 4 years. Our X-ray obser- 
vations can hence be used as a 'time machine' to directly 
measure the CSM density during the last stages of the 
stellar evolution. 

Fig. 2 illustrates the change in mass-loss rate as a func- 
tion of the stellar wind age. It can be seen that the 
mass-loss rate constantly declined from M = 4 x 10~ 4 
to 4 x 10~ 5 M Q yr _1 (u w /10 km s -1 ) just prior to the 
explosion. Integration of the mass lost by the progeni- 
tor along the shock front gives the CSM density at the 
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given interaction radius from the SN site. The CSM den- 
sity distribution is illustrated in Fig. 3. It is found that 
the CSM density profile is best described by a power law 
Pcsm °c r~ s with index s — 1.63 ± 0.02, assuming a spher- 
ically symmetric expansion. This is significantly flatter 
than expected for a constant mass-loss rate and constant 
wind velocity profile (r~ 2 ). The results confirm previous 
evidence for a rather flat CSM density profile during the 
first year (s ~ 1.5, van Dyk et al. 1994; 1.5 ^ s <J 1.7, 
FLC96; s = 1.661^5, Marcaide et al. 1997). 

The observed change in the pre-SN history could be due 
to variations in one of the following parameters: 

• Variations in velocity and/or temperature of the 
shocked CSM may have caused changes in the X-ray out- 
put of the shocked region. This model was proposed to 
explain the radio lightcurve of SN 1993J (Fransson & 
Bjornsson 1998). In fact, evidence for a ^ 20% decrease 
in the expansion velocity over the first ~ 4 years was re- 
ported (Marcaide et al. 1997; Bartel et al. 2000). This 
change in shock expansion velocity, however, would only 
lead to a < 10% increase in the mass-loss rate and cannot 
account for the observed order-of-magnitude change. Also, 
as the ROSAT soft X-ray band cooling function for opti- 
cally thin thermal plasma emission is not very sensitive 
to the plasma temperature (A(T) oc T 0,5 ; e.g. Raymond, 
Cox & Smith 1976), a change in temperature from 10 7 K 
to a few xlO 11 K would be required to explain the ob- 
served change in the mass-loss rate history. Such drastic 
temperature variations have neither been observed for SNe 
nor have they been put forward by models describing the 
CSM interaction. Therefore, variations in either the shock 
expansion velocity or in the temperature accounting for 
the inferred evolution can be ruled out. 

• A different scenario would be a non-spherically sym- 
metric geometry caused by a binary evolution of the pro- 
genitor. Observational evidence that the progenitor of 
SN 1993J was a stripped supergiant in a binary system 
has been presented based on the visual lightcurve, which 
cannot be described by a single, massive star alone (e.g. 
Podsiadlowski et al. 1993). In addition, "double-horned" 
emission line profiles indicate the presence of a flattened or 
disk- like expanding shell (Matheson et al. 2000a; 2000b). 
The X-ray lightcurve also shows some significant devia- 
tions from the long-term i~ - 27 behavior (cf. Fig. 1). 
Such deviations have also been observed for SN 1979C 
in the radio regime (Montes et al. 2000) and were dis- 
cussed in the context of a self-colliding binary stellar wind 
model leading to the formation of a multiple shell-like CSM 
profile (Schwarz & Pringle 1996). X-ray observations of 
SN 1979C, however, were not suited to address these ques- 
tions (Immler, Pietsch & Aschenbach 1998; Kaaret 2001). 
Similarly, a binary model for SN 1993J cannot be chal- 
lenged by the ROSAT data alone and no information about 
the long-term radio lightcurve of SN 1993J is available. 

• Alternatively, the observed evolution is caused by 
changes in the wind parameters. The X-ray lightcurve 
implies p C sm r~ s with s = 1.63, or, using the con- 
tinuity equation, M/v w cx r 2 ~ s oc r - 37 oc t^ 37 - The 
ratio M/v w therefore decreases with time approaching 
the date of the explosion, indicating either a decrease in 
the mass-loss rate, an increase in the wind speed or a 
combination of both. It is interesting to compare typi- 



cal wind parameters for red and blue supcrgiants, which 
are on the same evolutionary track for massive (<; 15M Q 
ZAMS) stars (e.g. Salasnich, Brcssan & Chioso 1999). 
Whereas the mass-loss rate at t w <~ 10 4 years (M = 
4 x 10~ 4 M yr~ 4 (?; w /T0 km s" 1 )) is typical for a red 
supergiant (e.g. M = 3 x 10~ 4 M yr" 1 for HD 179821, 
Jura, Velusamy & Werner 2001; Chin & Stothers 1990; 
Schaller et al. 1992), the observed change in ratio M/i> w 
by an order of magnitude at t w ^ 30 years could indicate a 
transition between the red and blue supergiant phase due 
to an increase in effective temperature of the star. Blue su- 
pergiants are known to have significantly lower mass-loss 
rates and higher wind velocities (M ~ 10~ 6 M Q yr -1 , 
v w ~ 500-1000 km s" 1 ; e.g. Kudritzki et al. 1999). 
This scenario for the evolution of the SN 1993J progen- 
itor would have some interesting similarities with that of 
SN 1987A, whose progenitor completely entered the blue 
supergiant phase after significant mass-transfer to a com- 
panion (Podsiadlowski et al. 1993). 

Whereas our X-ray data alone cannot give information 
on whether M or v w have effectively changed during the 
late stage of the progenitor's evolution, there is a clear dif- 
ference between these possibilities regarding the kinetical 
energy of the stellar wind: using p csm oc r - 163 ; the kineti- 
cal wind energy decreased by a factor of 13 for a changing 
M during the observed At w ~ 10 4 years. In the scenario 
of a changing v w only, an increase in kinetical wind energy 
by a factor of 10 22 is expected. Future stellar evolutionary 
models might give answers as to which of the two possi- 
bilities is more likely. 

4.2. X-ray emission from the shocked-heated SN ejecta 

Let us explore the scenario that the X-ray emission is 
dominated by the ejecta, heated by the reverse shock. A 
necessary condition for the ejecta accounting for the ob- 
served X-ray emission is that the absorption by the post- 
shock gas is low. In this model, an ejecta density struc- 
ture of p sn = po(t/to)~ 3 (v/v )~ n is expected (FLC96). 
The total X-ray output of both the reverse-shock heated 
ejecta and the shocked stellar wind in the forward shock 
is then given by L x oc (M/v w ) 2 x (t/t ) 3 ~ 2s in the case 
where n is large. In fact, using our inferred X-ray rate 
of decline of L x oc i~ - 27 is entirely consistent with the 
above L x oc t 3 ~ 2s oc i~ - 24 for s = 1.63. The assump- 
tion that n must be large is hence confirmed by our data 
and has also been concluded from radio and optical obser- 
vations of SN 1993J (n ~ 20-30; Shigeyama et al. 1994; 
Baron, Hauschildt & Branch 1994; Suzuki & Nomoto 1995; 
FLC96). It is important to note that if the density gradi- 
ent for the ejecta is large, X-rays from the reverse shock 
must be heavily absorbed. Whereas variations in tempera- 
ture of the ejecta can be ruled out to explain the observed 
X-ray lightcurve (cf. §4.1), a change in absorption from 
initially 10 22 cm~ 2 to the Galactic column of 4 x 10 20 cm~ 2 
could account for the observed decrease in source flux (cf. 
Fig. 5 in Wang, Immler & Pietsch 1999). The signature of 
high absorption, however, is absent in all ROSAT PSPC 
spectra of SN 1993J, which are consistent with the Galac- 
tic foreground absorption (Nh = 5.3 ± 1.7 x 10 20 cm~ 2 , 
Zimmcrmann et al. 1994, 1996). The most likely scenario 
is hence that the emission from the reverse shock is com- 
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pletely absorbed and that the observed soft X-rays are only 
due to the shocked CSM of the progenitor wind. 

Based on the modeling of the early X-ray lightcurve of 
SN 1993J, Fransson, Lundqvist & Chevalier (FLC96) con- 
cluded that the X-ray emission during the first months 
originates from the interaction of the SN shock with the 
circumstellar medium. As the reverse shock begins to pen- 
etrate the cool shell it is expected to contribute an increas- 
ing fraction to the total X-ray output. This predicted rise 
of the ROSAT band X-ray flux after ~ 100 days (as illus- 
trated in Figs. 8 and 10 in FLC96) due to the emerging 
ejecta is not observed with ROSAT. Instead, the overall 
long-term X-ray lightcurve (Fig. 1) is declining with a 
power-law close to the initial measurements, which were 
consistently explained as a result of the interaction of the 
SN shell with the ambient CSM. 

5. CONCLUSIONS 

We present the first detection of a prc-SN evolution in 
X-ray, based on long-term monitoring of SN 1993 J with 
ROSAT over a period of 5 years. The data are fully con- 
sistent with a description in the context of the SN shock 



interacting with the CSM blown off by the progenitor's 
stellar wind. From the X-ray rate of decline L x oc i~ - 27 
we infer a CSM profile p csm oc r~ 163 , which is signifi- 
cantly flatter than expected for a constant mass-loss rate 
and constant wind velocity profile {r~ 2 ). The observations 
cover an epoch of ~ 10 4 years in the progenitor's stellar 
wind history. During this period, the mass-loss rate of the 
progenitor has decreased constantly from M = 4 x 10~ 4 to 
4 x 10~ 5 M yr _1 (w w /10 km s _1 ) just prior to the explo- 
sion. The most likely explanation for this pre-SN evolution 
is either an increase in wind speed, a decrease in mass-loss 
rate or a combination of both, indicative that the progen- 
itor star was undergoing a transition from the red to the 
blue supergiant phase. The data demonstrate the scien- 
tific potential of long-term X-ray monitoring of SNe as an 
important diagnostical tool to probe the CSM interaction 
and the evolution of SN progenitors. 

This research made use of various online services and 
databases, e.g. ADS, HEASARC, NED, and the ROSAT 
data archive at MPE. The project is supported by NASA 
under the grants NAG 5-8999 and NAG5-9429. 
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2. — Mass-loss rate history of the SN 1993J progenitor. PSPC data arc marked by boxes, HRI data are indicated by crosses. 
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Fig. 3. — Circumstcllar density profile as a function of SN shell expansion radius. The solid line gives the best-fit CSM profile of p C sm oc 
to the PSPC (boxes) and HRI (crosses) data points. 



